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The production of new particles at a hadron collider like the LHC is always accompanied by 
QCD radiation attributable to the initial state (i.e. ISR). This tends to complicate analyses, so ISR 
is normally regarded as a nuisance. Nevertheless, we show that ISR can also be valuable, yielding 
information that can help in the discovery and interpretation of physics beyond the Standard Model. 
To access this information we will introduce new techniques designed to identify ISR jets on an event- 
by-event basis, a process we term ISR tagging. As a demonstration of their utility, we will apply 
these techniques to SUSY di-squark (di-gluino) production to show that they can be used to identify 
ISR jets in roughly 40% (15%) of the events, with a mistag rate of around 10% (15%). We then 
show that, through the application of a new method which we will introduce, knowledge of an ISR 
jet allows us to infer the squark (gluino) mass to within roughly 20% of its true value. 



I. INTRODUCTION 

Quarks and gluons are always splitting apart and re- 
combining [1 . When they are scattered at high energies 
this process is interrupted, and the result is that addi- 
tional quarks and gluons which could not be recombined 
end up in the final state. This radiation, attributable 
to the splitting of the incoming states, is termed initial 
state radiation (ISR). 

ISR often complicates analyses. For example, it can 
overlap with, and thus contaminate, jets formed from 
the decay of new particles (we will call these FSR jets). 
Furthermore, when ISR emissions yield additional inde- 
pendent jets (i.e. ISR jets), sorting out the combinatorics 
of an event can be even more difficult. These are not ir- 
resolvable difficulties, and indeed recently progress has 
been made toward removing sources of jet contamina- 
tion [2], mitigating combinatorial difficulties [3], defin- 
ing new observables less sensitive to contamination 1], 
and more accurately accounting for the physics of ISR in 
Monte Carlo simulations [5]. 

However, ISR - rather than always proving to be an 
obstacle - can actually be helpful in the study of physics 
beyond the SM. This has already been demonstrated in 
several recent studies. In Ref. [6], ISR was shown to make 
BSM events more prominent by giving pair-produced 
new-physics states something to recoil against, thus in- 
creasing both and Ht- Meanwhile, Ref. [7j took a 
different approach, calculating the influence of ISR on 
inclusive variables (which sum over everything measured 
in the detector) with the aim of inferring the mass scale 
probed in new-physics events. 

ISR's utility lies primarily in the fact that the physics 
of its production is (approximately) dependent only upon 



the mass scale of the event and the states it couples to 
- independently of any intermediate or final decay prod- 
ucts. So while the kinematic quantities measured of FSR 
jets in BSM processes depend upon some combination 
of all the masses involved in the process, the kinematics 
of ISR depend upon only the masses of the particles in- 
volved in the production process (i.e. not upon those pro- 
duced subsequently in a cascade/decay) and their cou- 
plings to initial state particles. 

Here we will make use of this information carried by 
ISR, with the main new feature of our approach being 
that rather than treating ISR in an inclusive way, as 
in Ref. [7 , we will seek to identify a particular jet as 
attributable to ISR. We will see that such jets can be 
identified on an event-by-event basis with a small mis- 
tagging probability. We will then show that this sort of 
technique for tagging ISR jets provides a powerful tool 
for understanding BSM physics that can be applied to 
many interesting processes. We focus in particular on 
a new kinematic method that can be used to evaluate 
squark and gluino masses. 

The outline of this paper is as follows. Sec. [IT] will 
discuss the sort of techniques one can use to tag ISR jets, 
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Sec. Ill will discuss one example of what we can learn 
about BSM physics through the study of ISR jets based 
solely on kinematics, and Sec. |TV] will apply these ideas to 
di-squark and di-gluinio production as a demonstration 
of their feasibility and utility. Finally, Sec. |V] contains 
our conclusions. 



II. TAGGING AN ISR JET 

Tagging an ISR jet requires identifying characteristics 
that distinguish it from the other jets in an event. Al- 
though ISR possesses some general traits that hold true 
regardless of the process at hand |8], tagging ISR jets 
solely based on these properties would be challenging, 
if not impossible. We instead focus on tagging ISR in a 
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particular class of interesting processes - the pair produc- 
tion of BSM particles, each of which decays into jets and 
an invisible particle (i.e. pp — >■ NjJ + 2x1 + where 
Nf = 2/4 for di-squark/di-gluino production). Although 
we will restrict ourself to these topologies, we expect that 
ISR jets are also identifiable in other cases and that sim- 
ilar techniques could be developed for more complicated 
processes. Nonetheless, this will serve as a proof of con- 
cept, that we employ later in Sec. |IV[ in what is already 
an important application of these ideas to BSM physics. 

It turns out that the modest assumption of pair pro- 
duction gives one a significant handle for identifying ISR 
jets. Suppose one expects to see Nf FSR jets in a BSM 
event. As long as there's no reason for these to be par- 
ticularly soft, one can assume that of the Nf + 1 hardest 
jets in the event, Nf are attributable to FSR and one 
to ISR. As the production process is symmetric, all of 
the properties governing the production of one FSR jet 
should hold for the others. Thus, of the Nf + 1 hardest 
jets in the event, we can identify the ISR jet as the one 
which is in some way distinguished from the others. 

The method we prescribe for accomplishing this is to 
consider the Nf + 1 hardest jets in an event and identify 
a candidate ISR jet (here labeled i) for which at least one 
of the following conditions is met [5] : 

1. The jet's pt is distinct (i.e. it is harder or softer 
than the others): 



max(pTi:,PTj) 
mm{pTi,PTj) 



> 2 V J ^ I 



(1) 



2. The jet is separated from the others in rapidity: 

\y,-y,\>1.5y J (2) 

3. The jet is distinguished by its rrii/pTi = ra- 
tio [TOl: 



max(Ai, Aj) 
min(Ai, Aj) 



> 1.5 V j 7^ i 



(3) 



If a jet (again labeled i) is selected by any of the above 
criteria it should then satisfy all of the following: 

• The selected jet must not be central: \yi\ > 1. 

• It must not be too close to the other jets, which are 
all implicitly FSR jets: 



IVi ~ Vjl > 0.5 V j ^ « 



(4) 



• These other jets must be reasonably close to each 
other in px'- 



PTj 
PTk 



1/2 

1 - a 



for PTj(k) = max(min){pTi|V / ^ i}, with p = 2(3) 
for Nf = 2(4), and where we have introduced the 
variable 



min (pti,^t) 



(6) 



(5) 



max(pTj,-^T) 
to relax this condition when the ISR is very hard. 

• Finally, the implicit FSR jets must be somewhat 
central: \yj\ < 2 V j 7^ z 

If any of the above conditions is not satisfied, the jet 
being considered is not tagged and other jets are checked 
to see if they pass any of the distinguishing criteria (Eqs. 

We note that it is surely possible to improve upon the 
technique presented above, and that the numerical values 
we presented have not been thoroughly optimized. Even 
so, we will see these criteria already work quite well, trig- 
gering on 40% (15%) of the events, for Nf = 2(4) topolo- 
gies, with a small 10% (15%) mistag rate. 



III. USES OF AN ISR JET 

Once an ISR jet has been identified in an event it can 
be used in multiple ways to shed light upon the under- 
lying physics that produced it. As the production of 
ISR is determined by the mass scale probed by the pro- 
cess, the identity of the partons in the initial state, and 
the relevant parton distribution functions (PDFs), the 
resulting ISR kinematical distributions will reflect all of 
these influences [Tl]. Here though, rather than focus on 
general properties of the aforementioned distributions, 
whose calculation would depend upon a careful treat- 
ment of QCD, we will instead present a simple new kine- 
matical technique useful in measuring ttibsM) the center 
of mass energy for the two heaviest BSM particles pro- 
duced in the symmetric processes we are considering. Be- 
cause hadron colliders tend to produce heavy states close 
to threshold, a measurement of tobsm is nearly equiva- 
lent to a measurement of the new-physics particle's mass: 

TOBSM = \/{Vql~g +Pq*/g)^ ~ 2m^/g. 

Other kinematic variables are also sensitive in some 
way to mBSM- Examples include M^g [T^], Mt^ [TB] . 
and their more advanced extensions |14) . Some recent 
works [15J have also made use of ISR to give their 
distributions additional structure. However, these tech- 
niques are in general sensitive to all of the masses in the 
decay chain, or only work for very specific processes (e.g. 
giuino stransverse mass |14j). 

Remarkably, by looking to ISR we can construct a new 
kinematic measure sensitive only to mesMj independent 
of any other assumptions on the spectrum. The basic 
idea behind this method stems from the observation that 
any BSM particles produced must be recoiling against 
ISR in the transverse plane. Boosting the FSR system 
back along the transverse plane to compensate for the 
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ISR jet's pt requires an assumption for the system's cen- 
ter of mass energy, and only when we have assumed the 
correct value will the boost function properly. Before 
proceeding, we note that while any BSM particles are, 
in fact, recoiling against all of the ISR particles (rather 
than only the leading jet), in practice the ISR jets assume 
a strong hierarchy and using only the leading jet to 
apply a boost will serve as a reasonable approximation. 

In detail, the method we prescribe to measure mesM 
using the ISR jet's kinematics is to: 
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Identify all of the visible FSR jets and boost them 
along the z direction so that the visible FSR is at 
rest in the z frame (i.e. the net Pz for FSR jets is 
zero). That is, each FSR four- vector (here labeled 
i) is shifted 



Pi 



+ Piz) 



(7) 



where (3 = —pz/E and 7 = — /J^, for pz 

and E the sum longitudinal momentum and en- 
ergy taken over all observable particles in the sys- 
tem. This boost is performed because, while ideally 
the system will be at rest in the z direction before 
boosting in the transverse plane (step two), this is 
a configuration we cannot achieve because of un- 
certainties introduced by missing energy. However, 
by applying the boost in Eq. ([t]) we approximate 
this condition. 

2. Boost the system along the direction transverse to 
the beam, parallel to the transverse momentum of 
the ISR jet, assuming some system mass M. This 
means that the projection of each FSR px vector 
along the ISR direction transforms as 




(8) 



where prp^ = px ■ p^^ is the projection of each px 
along the ISR pr direction. 

3. Measure the sum projection of the resulting 
boosted FSR along the ISR transverse direction, 
assigning the result a ±1 depending upon the sign: 



1-1 ifE,PT^<0 



(9) 



4. Finally, the average projection across many events 
is measured: (cr) = X^i^i ^i/-^ 

When (cr) is positive there is a net projection along the 
ISR axis, indicating the assumed mass is too small, while 
when it is negative, the assumed mass is too large. Ex- 
amples of the resulting distributions are shown in Fig. [l] 
for the case of di-squark and di-giuino production. 

Before proceeding, we call attention to two choices we 
made in the analysis that might be improved in a more 



TABLE I. ISR tagging efRciencies computed for different 
choices of spectra. The first four rows are for di-squark pro- 
duction, and the last four are for di-giuino production. 
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careful treatment. The first is in step one, where we 
boosted the FSR along the z direction to approximate the 
longitudinal rest frame. While this technique seems to 
operate reasonably well, it may be possible to better infer 
the z-boost using beam thrust techniques as suggested 
in Ref. [TB]. We further note that Eq. (|9| assigns each 
event an equal weight when computing (07, regardless of 
the measured imbalance. This choice was made because 
weighting events by their px imbalance (^^Pfi) tends to 
make (cr) sensitive to only a few outlier events. Perhaps 
a better measure exists, but we do not pursue it here. 



IV. EXAMPLE: DI-SQUARK & DI-GLUINO 
PRODUCTION 



We now apply the aforementioned techniques to the 
pair production of squarks and gluinos, letting q — >■ q+Xi 
and g qq + Xi (via an off-shell squark) . To perform 
this analysis we use Madgraph v4.4.51 [T7] to gener- 
ate 10^-event samples at matrix-element level, assum- 
ing a 14 TeV LHC, which are then showered in Pythia 
v6.422 [18 and matched using the MLM procedure 19 . 
Fully showered and hadronized events are then grouped 
into 0.1 X 0.1 cells (77, ip) cells between —5 < t] < 5, which 
are clustered in Fast jet v2.4.2 |20j using the anti-fcT- 
algorithm [21 . Our di-squark samples were clustered us- 
ing R — 0.7, while R — OA was used for the busier di- 
gluino events. Note that, to simplify matters, we have 
not accounted for the effects of multiple interactions or 
pileup. 

Table |l] shows the efficiencies found using the tagging 
procedure of Sec.|ll] Remarkably, we see that the tagging 
efficiency (i.e. the percent of events in which an ISR 
jet is tagged) and the mistag rate (the percent of events 
in which a jet that has been tagged as ISR was tagged 
incorrectly) are stable, even when comparing a standard 
SUSY spectrum with tolsp — 100 GeV to one in which 
the LSP is nearly degenerate with the supersymmetric 
particle that decayed into it. 

Fig. JTl shows the distribution of (cr) for the spectra in 
Table [irwhere we see that the kinematical technique in- 
troduced earlier works quite well: (tobsm) = 1-3 TeV 
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for rriq/mg = 500 GeV and (tobsm) = 2.5 TeV for 
mq/rrig = 1 TeV, values which are quite close to those 
measured by the point at which the FSR momentum have 
been boosted to have no preferred direction {{a) = 0). 

V. CONCLUSION 

While ISR is normally regarded as a nuisance, here we 
have seen that it can instead be useful, allowing for qual- 
itatively new measurements of BSM physics that would 
be difficult or impossible to otherwise access. In this pa- 
per, we have introduced a set of techniques for tagging 
the ISR created in the pair production of BSM states, 
each of which decays into jets and an invisible particle. 
Although the methods we introduced are specific to this 
sort of process, they can be readily extended to scenar- 
ios where BSM physics realizes a more complicated final 
state topology. 

We have applied these techniques to SUSY di-squark 
(di-gluino) production, where we saw they tagged ISR 
in roughly 40% (15%) of the events, with a mistag rate 
of around only 10% (15%). This allowed us to make a 



qualitatively new measurement, using a kinematic tech- 
nique we introduced, of the squark/gluino mass to within 
roughly 20% of its true value, without making any as- 
sumptions on the rest of the SUSY spectrum. 

Given the success of this relatively simple procedure, 
we expect that more elaborate ISR taggers (employing, 
for example, the information contained in jet substruc- 
ture) could be constructed for use in other topologies to 
produce new and complementary measurements of both 
SM physics and physics beyond the Standard Model. 



ACKNOWLEDGMENTS 

We would like to thank D. Feldman, E. Kuflic, I. Kim, 
M. Lisanti, G. Salam, M. Schwartz, I. Stewart, and J. 
Thaler for helpful discussions. DK is supported by a Si- 
mons postdoctoral fellowship and by an LHC-TI travel 
grant. LR is supported by NSF grant PHY-0556111. 
L.-T.W. is supported by the NSF under grant PHY- 
0756966, and by a DOE OJI award under grant DE- 
FG02-90ER40542 



, B126 



298 (1977)1 



[1] G. Altarelli and G. Parisi, ]NucL Phys.. 
[2] J. M. Butterworth, A. R. Davis on, M. Rubin^ and 
G. P. Salam, Phys.Rev.Lett., 100, 242001 (2008) 
arXiv:0802.2470 [hep-ph] ; S. D. Ellis, C. K. Ver- 
milion, and J. R. Walsh, Phys.Rev., D80, 051501 
(2009), arXiv:0903.5081 [hep-ph]; D81, 094023 (20TO)[ ' 
jarXiv:0912.0033 [hep-ph]; D. Kroh n, J. Thaler, and L.- 
T. Wang, JHEP, 1002, 084 (2010) , arXiv:0912. 1342 [hep- 



Ph]| 
[3] J. Alwall 



K. Hiramatsu, M. M. 

Phys.Rev. Lett., 103, 151802 
[hep-ph] 



Nojiri, and Y. Shimizu, 
(2009)1 |arXiv:0905.1201| 



M. 



Sakurai, (2010), 
K. Kong, K. T. 



|arXiv:1008.1813| 



[5] 



M. Nojiri and K. 

[hep-ph]; P. Konar, K. Kong, K. T. Matchev, and 
M. Park, (2010), arXiv: 1006.0653 [hep-ph] 
J. Alwall, S. Hoeche, F. Krauss, N. Lavesson, 
L. Lonnblad, et ai, Eur.Phys.J., C53, 473 (2008J^, 
|arXiv :0706.2569 [hep-phj^ 

and O. Mattelaer, (2010), l arXiv:1010.2263 

T. Plehn, D. Rainwater, and P. Z. Skands, 
B645, 217 (2007), arXiv:hep-ph/0510144" 
J. Alwall, S. de Visscher, and F. Maltoni, 



[hep- 



Phys.Lett 
r [hep-ph] 

[JHEP, 0902, 017 (2009), arXiv:0810.5350 [hep-ph 
[6] J. Alwall , M.-P. Le, M. Lisanti, and J. G. Wacker, 
B666, 34 (2008), arXiv:0803.0019 [hep-pji_ 
D79, 015005 (2009), |arXiv:0809.3264 [hep 



Phys. Lett.. 



IE 



Phys.Rev., 


A. Papaefstathiou and B. Webber, 
"(2009), arXiv:0903.2013 [hep-ph]; 



JHEP, 0906, 069 
1007, 018 (2010), 

]arXiv; 1004.4762 [hep-ph]." 
[8] E.g. it is not very correlated in rapidity with FSR jets. 
[9] It is poss ible, although rare, for more than one jet to pass 
Eqs.( 1|3 1. If this happens we rank the tagging criteria by 
the order in which they are listed, e.g. accepting a jet 



which passes Eq. ([T| over one which passes only Eq. ([3|, 
etc. If more than one jet pass the same tagging criteria, 
we take as our candidate for Eq. ([T]) the hardest jet, and 
for Eq. (|3| the jet with the largest A. If more than one 
jet pass Eq. ([2|, none are accepted. 
[10] We note that^ is sensitive to higher order emissions, in 
contrast to the criteria of Eqs.([l|2|. Further work along 
this line, particularly in employing jet substructure (see 
Ref. ^22} for a review), could potentially improve upon 
the methods presented herein. 
[11] We will pursue these ideas in a foUowup paper. 
[12] I. Hinchliffe, F . Paige, M. Shapiro, J. Soderqvist, 

and W. Yao, Ph ys.Rev., D55, 5520 (1997)] |arXiv:hep-| 

' ph/9610544 [hep-phj' 



"J Alwafi, A. Freitasj [ 



[13] C. Lester and D. Summers, |Phys.Lett., B463, 9 9 (1999)J 
arXiv:hep-ph/9906349 [hep-ph]; A. Barr, C. Lester, and 
_P. Stephens, J.Phys.G, G29, 2343 (2003), arXiv:hep- . 
ph /0304226 [hep-ph]) 
[14] C. Lester and A. Barr , 
[arXiv:0708.1 028 [hep-ph]) 

Kim, and C. B. Park, "Phys.Rev.Lett., 100 
(2008), arXiv:0709.0288 [hep-ph]; H.-C. Cheng, 
gelhardt, J. F. Gunion, Z. Han 
252001 (2008) 



JHEP, 0712, 102 (2007)] 



W. S. Cho, 



K. Choi, Y. G. 

171801| 
D. En- 

and B. McEhath, |1(30^ 

arXiv:0802.4290 [hep-ph] 

[15] P. Konar, K. Kong, K. T. M atchev, and M. Park, 
105, 051802 (2010), |arXiv :0910.3679| 
and K. M. Zurek, JHEPl 



Phys.Rev.Lett. 

hep-ph|| T. Cohen, E. Kuflik. 



1011, 008 (2010), arXiv:1003. 2204 [hep-ph] 
I. W. Stewart, F. J. Tackmann, and W. J. Waalewijn, 



[Phys.Rev.Lett. 
[hep-ph] 



TO'5, 092002 (20l0)] |arXiv:1004.2489| 



J. Alwall, P. 
quet, et al, 
[hep-ph]. 



Her- 



Demin, S. de Visscher, R. Frederix, M 

JHEP, 0709, 028 (2007)] |arXiv:0706.2"334 



5 



<0> 



0.08- 
0.06^ 
0.04 7 ' . 
0.02 ^ 

-0.02- +1 -.t 
-0.04^ 



..niLsp^lOOGeV 
-m,,,-450GeV 



1000 1400 1800 2200 2600 3000 

Mass guess [GeV] 







<a> 




mLS[.-100GeV 


A 0.1 
0.08 


1 

r 

- + 

- + , 








0.06 












0.04 


T 












0.02 




\ 












-0.02 




\ 
















+++ 


-0.04 












10 


00 14 


00 18 


00 22 


00 26 


00 30 


on 



Mass guess [GeV] 



AO.I- 
V 

0.08- 



0.06 7 
004 



1 



-0.02 - 
-0.04- 



-iii,,,MOOGeV 



1800 2200 2600 3000 
Mass guess [GeV] 



0.08 
0.06 
0.04 
0.02 


-0.02 
-0.04 7 

1000 1400 1S( 



,i,^100GeV 
„„^900GeV 



1 2200 2600 3000 
Mass guess [GeV] 



FIG. 1. The average sign of the FSR projection along the transverse ISR direction for, proceeding left to right, di-squark 
production using rriq = 500 GeV, rriq = 1 TeV, and then di-gluino production with mg — 500 GeV, rrig = 1 TeV, with the LSP 
mass indicated in the legends. The position at which the points intersect (a) = is what we would identify as ttibsm, i.e. it 
where the FSR momenta are balanced because the boost is 'correct'. We see that it is in general close to 2mq/g. Note that the 
errors indicated are just the statistical errors associated with our Monte Carlo sample sizes. 
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